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We show that Anderson locahzation in quasi-one dimensional conductors with ballistic electron 
dynamics, such as an array of ballistic chaotic cavities connected via ballistic contacts, can be 
understood in terms of classical electron trajectories only. At large length scales, an exponential 
proliferation of trajectories of nearly identical classical action generates an abundance of interference 
terms, which eventually leads to a suppression of transport coefficients. We quantitatively describe 
this mechanism in two different ways: the explicit description of transition probabilities in terms of 
interfering trajectories, and an hierarchical integration over fluctuations in the classical phase space 
of the array cavities. 
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I. INTRODUCTION 

The interplay of quantum phase coherence and re- 
peated random scattering is at the origin of many effects 
in mesoscopic physics.^ ^ These effects include weak local- 
ization and universal conductance fluctuations, both of 
which are small but fundamental corrections with respect 
to the conductance obtained from Drude-Boltzmann the- 
ory. They culminate in Anderson localization, the phe- 
nomenon that the resistance of a one or two-dimensional 
electron gas grows exponential ly w ith system size if the 
system size is sufficiently large 

Originally, Anderson localization and other mesoscopic 
effects were discovered in the context of disordered met- 
als, in which electrons scatter off impurities with a size 
comparable to their wavelength. Theoretically, quantum 
effects in disordered metals are described using the 'disor- 
der average' which deals with an ensemble of macroscop- 
ically equivalent but microscopically different impurity 
configurations. The presence of impurities is not essen- 
tial for the existence of quantum effects, however. The 
same effects, with the same statistical properties, have 
been found to appear if the electron motion is ballistic 
and chaotic, the only source of scattering of elect rons 
being specular reflection off the sample boundaries l^*^*''^ 
Besides being of theoretical interest for understanding 
the quantum properties of systems with chaotic classical 
dynamics,^ the case of ballistic electron motion is relevant 
experimentally for very clean artificially structured two- 
dimensional electron gases in semiconductor heterostruc- 
tures, such as quantum dots or antidot lattices. 

Unlike disordered metals, in which impurities scatter 
diffractively, electrons in ballistic conductors have a well- 
defined classical dynamics. Many quantum properties of 
ballistic conductors with chaotic classical dynamics can 
be understood in terms of classical mechanics by making 
use of well-chosen semiclassical methods.'^ A well-known 
example is the Gutzwiller trace formula, which relates 
the density of states to properties of periodic orbits of 



the classical dynamics.'^ For the conductance G and its 
quantum corrections, there is a variant of Gutzwiller 's 
formula, which expresses G as a double sum over classi- 
cal trajectories a and (3 connecting the source and drain 
contacts,!^ 

G - ^ y dpdg A„A^e^(•5"-5^)/^ (1) 

Upon entrance and exit, the trajectories a and (3 have 
the same transverse momenta p and respectively, but 
the position at which they enter or exit the sample may 
be different. Further, ^a,/3 and Sa,f3 are the so-called 
"stability amplitude" and the classical action of the tra- 
jectories. 

With Eq. ^ as a starting point, the conductance, in- 
cluding its quantum corrections, can be calculated solely 
from knowledge of classical trajectories. Quantum ef- 
fects have been linked to the existence of families of clas- 
sical trajectories that only differ near 'small-angle en- 
counters' of the trajectory with itself.^i^^® ^® This 
semiclassical approach has been very successful explain- 
ing quantum effects in the 'perturbative regime' in which 
quantum interferen ce provides a small co rrection to the 
classical mechanics. ESEMZEMmim^ ^eak localiza- 
tion and universal conductance fluctuations are examples 
of such perturbative quantum effects. Recently, Heusler 
et al. showed that it is also possible to understand quan- 
tum effects in terms of classical trajectories outside the 
perturbative regime. They considered the spectral form 
factor K{t) of a ballistic cavity, the Fourier transform 
of the two-point correlation function of the density of 
states. A direct semiclassical evaluation of K{t) using 
Gutzwiller's trace formula was known to be possible in 
the perturbative regime t < in only, where the Heisen- 
berg ti me = h//S., A being the cavity's mean level 
spacing.li^iSlloIIIIlllSl Motivated by the field-theoretical 
formulation of the problem, Heusler et al. used a differ- 
ent way to express K{t) in terms of periodic orbits that 
allowed them to calculate K{t) for all times.f^^ 
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In this communication, we consider Anderson localiza- 
tion for ballistic electron gases, and show that it, too, can 
be understood in terms of interference of classical trajec- 
tories, with Eq. ([T|) as a starting point. Examples of 
systems that may exhibit such 'ballistic Anderson local- 
ization' are antidot lattices or arrays of chaotic cavities. 
Although it is generally accepted that Anderson localiza- 
tion exists irrespective of the details of the microscopic 
electronic dynamics,!^ the similarity of the phenomena 
for ballistic and disordered electron systems should not 
obscure the vastly different starting points of the the- 
ories for the two cases. This difference not only per- 
tains to the microscopic dynamics (classical-deterministic 
vs. quantum-probabilistic), but also to the statistical as- 
sumptions of the theory. Unlike theories of quantum 
transport in disordered metals, semiclassical theories of 
ballistic conductors are intended to describe one specific 
system.fSl Fluctuations appear solely from variations of 
the Fermi energy; no changes in the classical dynamics 
are invoked. 

For disordered metals, Anderson localization is most 
prominent for a quasi-one dimensional geometry, with 
sample length L much larger than the sample width W . 
For quasi-one dimensional samples, a full theory of trans- 
port in the localized regime was developed by Dorokhov^ 
and Mello, Pereyra, and Kumar,'22] using a stochastic 
approach, and by Efetov and Larkin,'221MJ using a field- 
theoretic approach. Our theory of ballistic Anderson lo- 
calization closely follows these approaches. The fact that 
the semiclassical theory for ballistic electrons follows the 
corresponding quantum mechanical theory for disordered 
metals is not special to the present case. It is also typ- 
ical of the trajectory-based semiclassical theories in the 
perturbative regime, which have a structure that resem- 
bles the diagrammatic perturb ation th eory of quantum 
corrections in disordered metals.^ ^ * ^^ * "'^^ * 

In addition to the ballistic electron gases considered 
here, Anderson localization also occurs in certain dy- 
namic systems with a periodic time-dependent Hamil- 
tonian, such as the kicked rotor.^I Classically, the dy- 
namics of the kicked rotor is chaotic, with a momentum 
that changes diffusively under the influence of the peri- 
odic kicks. Quantum mechanically, it exhibits 'dynamic 
localization', Anderson localization in momentum space 
rather than in real space.l^ The phenomenology of dy- 
namic localization is equal to that of Anderson localiza- 
tion in disordered met als. T his is confirmed by extensive 
numerical simulations j^^^ as well as a field theoretical 
analysis of the problem.l^ 

The trajectory-based theory of ballistic Anderson lo- 
calization that we report here is constructed for a spe- 
cific system, an array of chaotic cavities. A schematic 
drawing of such an array is shown in Fig. [T] We use 
this model system because the trajectory-based theory 
of transport through a single chaotic ca vity is well estab- 
lished in the Hterature.'SElIiMlIIIIlZIll Arrays of cavities 
have been used as a starting point for a field-theoretic de- 
scription of Anderson localization using random matrix 




FIG. 1: (Color online) Schematic drawing of an array of n 
chaotic cavities. The cavities are connected via ballistic con- 
tacts with dimensionless conductance Qc- The semiclassical 
theory requires the limit c/c — > oo at a fixed ratio n/^c- 



theoryj^Slini but not with Dorokhov's method. In Sec. |TT] 
we show how Dorokhov's theory can be adapted to this 
system if the cavities are disordered and random ma- 
trix theory can be used to describe transport through 
a single cavity. In Sec. Ill we summarize the basic el- 



ements of the trajectory-based semiclassical formalism. 
This formalism is used to construct the trajectory-based 
semiclassical theory of ballistic localization in Sec. |IV[ 
In section |V] we approach the localization phenomenon 
from a different perspective. We describe the array in 
terms of a nonlinear cr- model whose perturbative ('di- 
agrammatic') evaluation generates structures of paired 
Feynman amplitudes similar to those appearing in the 
native semiclassical approach. Alternatively, the dynam- 
ical structure of the theory can be analyzed to identify 
and successively eliminate hierarchies of different types of 
dynamical fluctuations in the system. In this way we ob- 
tain an effective low energy theory which turns out to be 
equivalent to the nonlinear cr-model of diffusive quantum 
wires; the latter model is known to predict exponential 
localization at large length scales. We conclude in Sec. 

m 



II. ARRAY OF DISORDERED CAVITIES 

We now describe how Dorokhov's approach can be used 
to describe Anderson localization in an array of disor- 
dered cavities or quantum dots. We take the disorder in 
each cavity to be weak (cavity size much smaller than 
the localization length) , so that transport through a sin- 
gle cavity is described by random matrix theory.^^ 

A schematic drawing of the array of cavities is shown 
in Fig. [T] The cavities are connected via ballistic con- 
tacts with dimensionless conductance gc- Since we want 
to compare with a semiclassical theory for an array of 
ballistic cavities, we require 3> 1. Localization takes 
place if the conductance of the array is of order unity. 
This condition is met if the number n of cavities in the 
array is comparable to gc. Hence, in the calculations be- 
low we take the limit ^c ^ oo while keeping the ratio 
n/ fixed. The same limit is take n in the field theo- 
retical description of localization, '^^^--^IinilinEl] -^j-^gj-g ]^ jg 

known as the "thick wire limit" . 

The transport properties of the array of cavities are 
encoded in its scattering matrix S'^{n). (The superscript 
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"q" is used to distinguish the quantum mechanical scat- 
tering matrix from its semiclassical counterpart to be in- 
troduced in Sec. Ill ) The matrix indices of S'^{n) rep- 
resent the two contacts « = 1,2 at the far left and right 
of the array and the transverse modes in each contact, 
\p±^i\ = ztTrhrii/W, Ui — 1,2, . . . , N, where N = gc is 
the number of channels in contact i, i = 1,2. The ma- 
trix S"^(n) is a random quantity because it depends on 
the Fermi energy and on the precise disorder configura- 
tion in each cavity. Followin g the approach of Dorokhov, 
Mello, Pereyra, and Kumar,'2S129| -^q consider the hermi- 
tian matrix 



(2) 



and calculate its statistical distribution by expressing 
T''(n) in terms of T'^(n — 1) and proceeding recursively. 
The matrix T'^ is related to the dimensionless conduc- 
tance g{n) of the array through the Landauer formula. 



g{n) ^ trr'i(n). 



(3) 



Taking the scattering matrix of each individual cavity 
from the circular ensemble of random matrix theory,Hil 
one then finds that the recursion relation for T'^ takes 
the form 



^ T'^{n) ~T\n-l) 

= -—T'^(n-\)tYT'^(n-\) 
9c 

+ X^{n) + 0{g-^/'), 



(4) 



where the hermitian matrix A" is a random (noise) term 
with a Gaussian distribution. 



(5) 
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where T'^ = T'^{1 — T'^), /3 = 1 or 2 in the presence or 
absence of time-reversal symmetry, respectively, and 



Qq _ cq cqt c-q 



(6) 



The averaging brackets (...)„ denote an average with re- 
spect to the disorder configuration in the nth cavity only. 

In the limit gc — » oo while keeping n/gc fixed, 
the stochastic recursion relation Q can be mapped to 
the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation, 
which is a stochastic differential equation for the eigenval- 
^qmrnEnm^he solution of the DMPK equation 
is known,'Sl which completes the theory of localization 
for an array of disordered cavities. 

Alternatively, the stochastic recursion relation (4| can 
be used to generate a coupled set of recursion relations 
for the disorder averages of the moments of T(n), 



\m— 1 



m—1 I \m— 1 / 

- - E E - ^^-^+1)) n ) 



k=l J=l 

n k — 1 



ri = l 



E E ( (^"^+»^ - n ) + ^^9c^)^ 



m = l 



(7) 



with 

T„ = tr(rq)". (8) 

(The argument n — 1 is suppressed on the right hand 
side of the second equahty.) The average in Eq. ^ is 
the full disorder average, applied to all cavities in the ar- 
ray. Taking the limit — > oo at fixed n/gc, Eq. ^ is 



mapped to a coupled set of differential equations for the 
moments of T'^ which is identical to the corresponding 
set of differential equations for a disordered wire .1^^^ A 
subset of these equations can be resummed into a par- 
tial differential equation for the generating function Ff^, 
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witlPEH 

fi = (detH 



2+ (cosCgg) - i)Tq 
2+(cosh(0i±6i2)-l)Tq 



l/2\ 



F, = < det 



2 + (cos(g3)-l)rq 
2+ (cosh(6'i) - l)Tq 



(9) 



As shown in Refs. 140 146 1 the resulting theory of local- 
ization in quasi-one dimension is formally equivalent to 
that o btained fr om the one-dimensional nonlinear sigma 
model !3" l 3i | 39 l 4 ll 



III. SEMICLASSICAL FORMALISM 

The central object in the trajectory-based semiclassi- 
cal theory of localization in an array of ballistic chaotic 
cavities is a semiclassical representation of the scattering 
matrix 5^. In the semiclassical representation, the dis- 
crete transverse momenta become continuous variables, 
so that the scattering matrix S'^ becomes a 'scattering 
kernel' Sij{p',p). Following standard semiclassical ap- 
proximations, this scattering kernel is then represented 
as a sum over clas sical trajectories a connecting contact 
j to contact 



Stj{p',p) 



1 



(10) 



such that the transverse momentum of a upon entrance 
and exit equals p and p' , respectively. Further, Sa is the 
classical action of a, and the stability amplitude. 



A- 



dq 



-1/2 



(11) 



where q is the transverse position upon entrance into the 
sample. Maslov indices and other additional phase shifts 
are included in Sa- Because the transverse modes in the 
quantum mechanical formulation are linked to the abso- 
lute value of the transverse momentum, not to the trans- 
verse momentum itself, the semiclassical counterpart of 
the products SijS\j and S^/^Sji consist of two contribu- 
tions: one in which transverse momenta in contact j are 
equal and one in which the transverse momenta are op- 
posite, 

^ {p , p")Slj {p, ±p") , 



£^sli±p",p')s^,{p",p), 

(12) 



Together, Eqs. (10 1 and (12 1 specify how products of the 



conjugate are expressed in terms of classical trajectories. 
[The "— " terms in the summations were omitted from 
the semiclassical expression for the conductance, Eq. ^ 
above.] 

For a theory of Anderson localization, we are interested 
in the trace of a product of alternating factors S and 
or in the product of s uch traces. Using the semiclassi- 
cal representation (10), a polynomial function F{S,S^) 
that involves the alternating product of n factors S and 
n factors is written as a summation over 2n classi- 
cal trajectories ai, . . . , a„ and . . . , one trajectory 
for each factor S or S"^ , respectively. Each configura- 
tion of classical trajectories is weighed by a phase factor 
exjp(iAS/h) with 



AS = Y,Sa,~Y.^0 



(13) 



quantum-mechanical scattering matrix and its hermitian 



Building on work by Richter and Sieber,liSliSliIl Haake 
and coworkers have identified a hierarchy of families of 
classical trajectories ai,. . . ,/3„ that contribute to the av- 
erage (F),^^ ^■^ where the average is taken with re- 
spect to variations of the Fermi energy while keeping the 
classical dynamics (i.e., the shape of the cavities) fixed. 
Their identification is based on the recognition that fam- 
ilies of trajectories contribute to (F) only if their total 
action difference AS is of order h systematically, which 
happens only if the trajectories ai are piecewise and pair- 
wise identical to the trajectories Pj, i, j = 1, . . . , n, up to 
classical phase space distances of order fi^/'^\lSllS\ Trajec- 
tories that are separated by larger phase space distances 
have typical action differences AS that are parametri- 
cally larger than h, so that their contribution vanishes 
upon taking the average. 

The simplest choice for a family of trajectories for 
which AS is of order h systematically is if each trajec- 
tory ai equals another trajectory Pj for the full length 
of the trajectory. Calculating (F) from the contribu- 
tion from such families of trajectories only is known as 
the "diagonal approximation" ) i4 | 37 | 38 | Nontrivial families 
of trajectories emerge from the possibility of small-angle 
encounters between trajectories, at which more than two 
trajectories are within a phase space distance ~ h}^'^\^At 
such small-angle encounters, the pairing between the ai 
and (3j can be changed — so that now trajectories need to 
be piecewise identical only. The duration of a small-angle 
encounter is the "Ehrenfest time" te — X^^ln{ppl/h), 
where A is the Lyapunov exponent of the classical dy- 
namics, pp the Fermi momentum, and / a characteris- 
tic length scale of the classical dynamics. The funda- 
mental action integra ls correspo nding to each small-angle 
encounter are known, J^IETMI g^^^ ^t^^ resulting theory 
takes the form of simple combinatorial rules with which 
any product of traces of products the scattering matrix 
and its hermitian conjugate can be calculated to arbi- 
trary order in H from the semiclassical representation of 
S, provided that the Ehrenfest time be much smaller 
than the sample's mean dwell time 7-j-, .[iMMSiga, (xhe 
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case of finite te/td is considerably more complicated, 
see, e.g., Refs. I24)25|49)50) but not relevant for a semi- 
classical theory of Anderson localization.) 

In the remainder of this text, we refer to a calculation 
of the energy-average {F) using contributions from fam- 
ilies of trajectories thus constructed as the "trajectory- 
based semiclassical formalism" . Although there is no for- 
mal proof that this formalism is exact, i.e., that there are 
no other contributions to {F) than from families of piece- 
wise paired classical trajectories, the formalism satisfies 
all known conservation rules and calculations based on 
trajectory-based semiclassics have been found to agree 
with fully quan tum mechanical calculations whenever 
applicable! ^^ * ^" * The present calculation can be viewed 
as another demonstration of the validity of trajectory- 
based semiclassics, by showing that the same formalism 
can serve as the starting point of a theory of localization. 

While we do not need the detailed results of the 
trajectory-based semiclassical formalism, there are two 
properties of ensemble averages calculated using that for- 
malism that are particularly relevant for our calculations 
below: 

(i) All averages are compatible with the condition of 



unitarityj 



[181 



3 3 

= SikS{p-p'). 



(14) 



(ii) For a product S^S\j or Sj^Sj)., the trajectory a of 
the semiclassical representation for S and the trajectory 
f3 of the semiclassical representation of S"^ at contact j 
satisfy 



Pl3, ka-q/sl^fi/PF 



(15) 



In particular, this implies that there is no contribution 
from the second term in Eq. (12 1 for a product of two 
scattering kernels.'^^Bsl 



IV. ARRAY OF BALLISTIC CAVITIES 

In the perturbative regime, the semiclassical theories 
of quantum corrections to transport and to the density of 
states closely followed the corresponding theory for dis- 
ordered metals. The semiclassical formalism described in 
the previous section was instrumental in formalizing the 
relation between the two types of theories. Motivated 
by this correspondence, we now look for the possibility 
to adapt the theory of localization in an array of disor- 
dered cavities [Sec. |ll] to the case of an array of ballistic 
cavities. 

Thus, paraphrasing the arguments of Sec. [ll] the goal 
of our calculation will be to find the full probability dis- 
tribution of the function 



for an array of n cavities. As shown in Sec.|lT] there are 
two ways in which this can be accomplished: 

1. Using a stochastic approach, in which one considers 
the stochastic evolution of the function T{n;p',p) 
as a function of n, or through 

2. the construction of a set of recursion relations for 
all moments of T(n). 

In both cases, the resulting theory is formally identical to 
the known theories of localization in quasi-one dimension. 

Although technically simpler, the stochastic approach 
is at odds with the goals of a semiclassical theory for 
ballistic localization: The goal of a theory of localization 
in an array of ballistic cavities is to describe an array of 
cavities with a fixed shape, using variations of the Fermi 
energy as the only source of statistical fiuctuations. Since 
the Fermi energy is set globally, for all cavities at the 
same time, quantum corrections for different cavities are 
not independent, and a stochastic approach is ruled out 
a priori. A stochastic approach is possible, however, if 
one relaxes the goals of the theory, allowing for small 
variations of the shape of each cavity, or for variations of 
a "gate voltage" that sets the Fermi energy of each cavity 
independently. 

Below we first describe the stochastic approach. In 
Sec. IIV Al we consider the case of broken time-reversal 
symmetry, which is technically simpler. The discussion 
of localization in the presence of time-reversal symmetry 



is given in Sec. IV B In Sec. IV C we discuss how a hier- 
archy of recursion relations for the moments of T can be 
constructed, where the average is taken with respect to 
variations of the Fermi energy of the entire array only. 



A. Stochastic approach 

The stochastic approach deals with (statistical) prop- 
erties of the function T(n) before averaging. Although 
the properties (i) and (i i) o f the trajectory-based semi- 
classical formalism [Sec. |III| are satisfied for the average 
of any product of traces of products alternating factors 
S and to arbitrary order in h, they have not been 
shown to follow from the semiclassical scattering matrix 
( 10 1 before averaging. However, since our goal is a sta- 

the final statements 



tistical theory of the transport 
of the theory will refer to averaged quantities only — 
we will accept these two properties on the level of the 
sample-specific semiclassical scattering kernel S{n;p',p) 
in the arguments that follow below. 

Starting from Eq. (10), the kernel T is expressed as 



a double sum over classical trajectories a and /3 that 
connect the entrance and exit contacts [Fig.|2], 



T{p\p)^ 



An A 



^ 27Th 



(17) 



T{n;p',p) = [Si2Sl2]{p',p) 



(16) 



Here p and q are the transverse momenta of a and f3 
upon entrance. The two trajectories have equal trans- 
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/3 



FIG. 2: (Color online) Example of a pair of trajectories a 
and j3 contributing to T{p',p) = [Si2Sl2]{p' ,p). Upon enter- 
ing the sample, the two trajectories have different transverse 
momenta p and q, and may enter at different transverse po- 
sitions. Upon exiting the sample, a and /3 have the same 
transverse momentum, and exit at transverse positions a dis- 
tance ^ h/pF apart. 



verse momenta upon exit, and exit at positions a quan- 
tum uncertainty ~ h/pp apart — see property (ii) above. 
Below, we express the difference ST = T{n) — T{n — 1) 
in terms of classical trajectories. We first calculate the 
average {ST)n, where the average is taken with respect 
to variations of the Fermi energy of the nth cavity (or 
of its shape), while keeping the Fermi energy and shape 
of the other cavities fixed. After that, we calculate the 
variance of dT and the higher cumulants. 

Average of ST. When calculating the average {ST)n = 
{T{n))n — T{n — 1) to leading order in 1/gc in the ab- 
sence of time-reversal symmetry, it will be sufficient to 
calculate (T(n))„ in the diagonal approximation by con- 
sidering trajectories a and /? that are "paired" in the 
nth cavity. Note, however, that trajectories need not be 
paired in the first n — 1 cavities, because no average is 
taken there. Each trajectory is classified by the number 
of times rria, mp that it enters the nth cavity from the 
(n — l)th cavity. Hence, for each trajectory there are 
TOq, and segments in the nth cavity, which we label 
as «!, . . . , and /3i, . . . , /3m^- Since trajectories are 
paired in the nth cavity, TOq. = mp = m. Examples of 
trajectory pairs a and [3 with m = 1, 2, and 3 are shown 
in Fig. [3] Since trajectories need to be paired upon exit, 
am has to be paired with While there are (ni — 1)! 
ways in which the remaining segments can be paired, we 
now show that the "diagonal pairing", aj paired with 
(3j, j = 1, . . . ,m, gives the main contribution to {ST)n, 
whereas all other pairings give contributions a factor 1 /gc 
smaller. 

Cutting the trajectories a and f3 at the contact be- 
tween the (n — l)st and nth cavity also separates the 
part of each trajectory that resides in the first n — 1 
cavities into m segments. The first segment of each tra- 
jectory connects the entrance contact to the exit of the 
(n— l)st cavity; All other segments connect the exit con- 
tact of the (n — l)st cavity to itself. Since the electron 
dynamics in the nth cavity is fully ergodic, the positions 
and transverse momenta with which these segments cross 
the interface between the (n — l)st and nth cavity are 
fully random, without correlations between different seg- 
ments. [Correlations can be ruled out down to quantum 
phase space distances ^ h/{ppLy/'^ because the dwell 





FIG. 3: (Color online) Examples of trajectories a and /3 
contributing to {T{n))n, where an average is taken over the 
Fermi energy of the nth cavity only. The four panels show 
trajectories contributing for m = 1 (top left), m = 2 (top 
right), m = 3 with diagonal pairing (bottom left), and m = 3 
with non-diagonal pairing (bottom right). The number m 
counts the number of times a and /? enter the nth cavity 
from the (n — l)st cavity. The three dashed lines in each 
panel represent the entrance contact, the contact between the 
(n— l)st cavity and the nth cavity, and the exit contact (from 
left to right). In the nth cavity, the trajectories a and (3 are 
piecewise equal up to quantum uncertainties. 



time td 3> te.] Hence, these m segments can be inter- 
preted as the semiclassical representation of a product of 



5i2, SI2, and m 



1 factors S22 and S'lj, before taking 



any average. Note that, although trajectories a,/? that 
are "paired" in the nth cavity have a phase space distance 
of order h/{ppLY^^ or smaller when they enter/exit the 
nth cavity from/to the (n — l)st cavity, such trajectory 
pairs are sufficient for the semiclassical calculation of the 



complete kernels T because of property (ii) of Sec. HI 

For the diagonal pairing, the m segments of a and 
(3 that reside in the first n — 1 cavities generate a par- 
ticularly simple product of factors S and S\ T{n — 
I)[tr7^(n - where 

7^(n;p,p') = [S22{n)^S22{n)]{p,p') (18) 

is the reflection coefficient for the first n — 1 cavities, seen 
from the exit contact of the (n — l)st cavity, and 



PF 



(19) 



tr7?.(?i) = / dpTZ{n]p,p). 

Using the trajectory-based semiclassical formalism to 
evalu ate th e diagonal trajectory sums in the nth 
cavity,'i^E^^'^ one then easily finds that the diago- 
nal pairing gives 

- E:^;^^(^-l)[tr7^(n-l)]'"-^ 

m=l ^ 5c 



_ T{n - l)gc 

~ 25c-tr7^(n- I)' 

Using unitarity, one has 

tr7?.(n) = 5c ~ tr T(n), 



(20) 
(21) 
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so that Eq. ( 20 ) can be rewritten as 



(<5r). 



1 



— r(n- l)trr(n- 1), 



(22) 



up to corrections of order l/g^- This is precisely the 
semiclassical equivalent of the recursion relation of Eq. 
Q, averaged over the Fermi energy or shape of the nth 
cavity. 

It remains to show that non-diagonal pairings, in which 
a segment aj is not paired with Pj give a contribution 
to {ST)n that can be neglected in the limit of large gc- 
Hereto we first consider the case m = 3, for which the 
only possible non-diagonal pairing is ai <-> /32, 0^2 ^ 
(3i [Fig. [3] bottom right]. For this pairing, the three 
segments of the trajectories that reside in the first n — 
1 cavities generate the semiclassical representation of a 

product of six scattering matrices, 'S'i2<S'|2'S'22'S'22'^22'5'i2- 
From unitarity, one has 



*5'l2 '^'22 "^22 "^22 '5*22 'S'12 



+ T{n-1)\ (23) 



For comparison, the diagonal pairing for m — 3 generates 
T{n — 1) [trTZ{n — 1)]^, which is a factor ~ larger 
because tr7?,(n— 1) ^ g^. Using the ergodic dynamics in 
the nth cavity, the non-diagonal pairing of segments in 
the nth cavity for rn = 3 gives a contribution to ((5T)„ 
equal to 



(ST) 



(3,non diag) 



^[T(n-l)-2T(n-l) 
+ T(n-1)3]. 



This is a factor ~ ^/gc smaller than the leading contri- 



bution (|22]) to {ST)„ . 

The same arguments can be used for m > 3: Non- 
diagonal pairings come at the expense of at least two 
factors tr7?.(n — 1) and, hence, lead to contributions to 
{5T)n that are at least a factor ~ l/gc smaller than the 
contribution from diagonal pairing. These arguments can 
also be used to show that contributions to (T(n))„ that 
involve small-angle intersections of the trajectories in the 
nth cavity are a factor 1/ gc smaller than the leading con- 
tribution considered above. 

Fluctuations of 5T . The fluctuations of 5T are de- 
scribed by the covariance {ST{pi,p2)6T{p[,p'2))n- We 
calculate {5T{pi,p2)ST{p[,p2))n from the identity 

{STsr),, = (5(rT'))„-r(n-i)(5r')„ 

-T'(n-l)(,5r)„, (25) 

where we used the shorthand notation T(n) = 
T{n\pi,p2), T'{n) = r(n;pi,p^), and 5{TT') = 
T{n)T'{n) - T{n - l)T'(n - 1). The average {5T)n 
is given by Eq. ( [22| above, so it remains to calculate 
{S{TT')). The product T{n)T'{n) is represented as a 
sum over four classical trajectories, a, (3, a', and (3' . As 
before, we introduce the numbers nia, mp, ma' and m^/ 





FIG. 4: (Color online) Examples ol four trajectories a, (3, a' , 
and /3' contributing to {T{n;p' ,p)T{n;p" ,p'))n for an array 
of chaotic cavities, where the average is taken with respect 
to the Fermi energy of the nth cavity only. The examples 
shown here have a non-diagonal pairing with ma = m^/ — 2, 
ttIq/ = 7710 = 1 (left panel) and ma = m^/ = rria' ~ m/s — 2 
(right panel). In the nth cavity, the trajectories a and a' are 
piecewise equal to (3 and l3' , up to quantum uncertainties. 



that indicate how often each trajectory enters the nth 
cavity. Since trajectories always enter or exit the nth 
cavity in pairs, one has TOq. -I- TOq' = -t- m^/ . 

Unlike the average ((5T)„, for which the only contri- 
bution came from the diagonal approximation in the nth 
cavity with diagonal pairing of the segments aj and Pj , 
the average of the second moment (((5TT'))„ has contri- 
butions from both non-diagonal pairings in the diagonal 
approximation and from trajectories beyond the diago- 
nal approximation, which have small-angle encounters in 
the nth cavity. We first consider the diagonal approxi- 
mation, for which each segment a; or a[ is paired with 
another segment (3j or (3j. For the diagonal pairing of 



segments, a. 



j = 1, . . . , nia' ~ 'm.pi , we find 



mp and a' ^ /?', 



(24) {5{TT'))f-^^T{n-l){5T')r,+T'{n-l){5T)n. (26) 



This contribution to {5T5T')n precisely cancels the sec- 
ond and third terms in Eq. (251. Hence {STST')n must 



be from non-diagonal pairings of the trajectory segments 
within the diagonal approximation, or from trajectory 
configurations beyond the diagonal approximation in the 
nth cavity. Since the latter class of trajectories have 
small- angle encounters, we write 



{ST6T'),, = {6{Tr)) 



non ciiae; 



(27) 



The leading non-diagonal pairing appears if one pairs 
the first u segments of a with the first u segments of (3' , 
and the last TOq — u segments of a with the last nia — u 
segments of f3, as well as the first v segments of a' with 
the first V segments of /?, and the last TOq' — v segments 
of a' with the last rria' — v segments of /?', where u = 
I, . . . , rria — 1 and v — 1, . . . , m'^ — 1 are integers. Two 
examples, with rria — 2, rriai = 1, u = 1, and v — 0, and 
with TUa = = 2 and u = w = 1, are shown in Fig. [4] 
One then finds 



(STST) 



non ciiae 



= —[^{Pl,P2)'T{pi,P2) 

9c 

+ T{pup'2)T{p\,P2) (28) 
+ T{p^,p'2)T{p'„p2)]+0{\lgl), 
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FIG. 5: (Color online) Two examples of a set of four trajecto- 
ries contributing to {TT')n that have a small-angle encounter 
in the nth cavity. For the left panel the encounter is in the 
interior of the nth cavity; For the right panel it touches the 
exit contact. 



where T(pi,p2) = T{n-l;pi,p2) and T{pi,p2) T{n- 
l;pi,P2), with 



T{n) = trS'i2S'22S'22'S'j2 
= T(n)-r(nf. 



Other pairings 
smaller. 



give contributions of order 1 / or 



(29) 



The second contribution to the fluctuations of 5T 
comes from trajectories with a small-angle encounter in 
the nth cavity. Since we only need fluctuations of 5T to 
leading order in 1/gc, it is sufficient to consider trajec- 
tories with one small-angle encounter only. Taking the 
small-angle encounter to be between the segments a„, 
(3'^, a'^, and /?(,, with 1 < u,v < m,m', and pairing the 
first It — 1 segments of a with the first u — 1 segments of 
/3' , and the last m — u segments of a with the last m — u 
segments of /3, as well as the first v — 1 segments of a' 
with the first v — 1 segments of /?, and the last m' — v 
segments of a' with the last m! — v segments of /?', we 
find 



oo oo oo 



{5T5T)T - -r(pi,p^)r(K,p2) E E E 



u.v—1 m—u 



1— / o^^^+^^' + ia^^"^'""*" 



ra —V 



^ ^ (tr 7^)"^^"^'^^ 



T{p,,p',)T{p[,p2)Y. E 2™+™'„«'-i 



--Tip,,p'2)np[,p2). (30) 
5c 



where the first line comes from encounters in the interior 
of the nth cavity and the secon d line comes from encoun- 
ters that touch the exit contact .I^HSil Examples of the two 
terms for m = m' — 1 are shown in Fig. [s] [We do not 
need to consider encounters that touch the contact be- 
tween the (n — l)st cavity and the nth cavity, because 
the contribution from such encounters is included in the 
(products) of the kernel T of the first n — 1 cavities.] 
Combining everything we have 

(<5T(pi,p2)<5r(p'l,P2)>n = -[.F(K,P2)T(pi,p^) 

9c 

+ T{pup'2)'T{p[,P2)] 

+ 0{\lgl). (31) 

Equation (31 1 is the semiclassical equivalent of Eq. ([s]). 

Higher cumulants of 5T can be calculated in the same 
way. For the fcth cumulant, one finds that only pairings 
that involve trajectories out of all k factors T contribute. 
Each additional factor T involved in the pairing scheme 
contributes an additional factor l/^c, which is why all 
cumulants with k > 2 are of sub-leading order in l/gc- 

Together, Eqs. (22) and (pTj) form the semiclassical 



equivalent of the stochastic recursion relation ([4| used 
in the fully quantum mechanical theory of localization in 
disordered quasi-one dimensional wires. In the limit g^ — > 
oo at fixed n/gc, Eqs. (22 1 and (311 can be mapped to 



stochastic differential equation for the eigenvalues of T, 
which will be formally equivalent to the DMPK equation. 
The DMPK equation, in turn, provides a full descrip tion 
of localization in quasi-one dimensional wires.^^SlSlIIU 



B. Presence of time- reversal symmetry 

In the presence of time-reversal symmetry both the av- 
erage {5T)n and the covariance {5T5T')n are different. 
In both cases the difference appears because one can pair 
time-reversed trajectories when taking the ensemble av- 
erage in the nth cavity. 

There are two additional contributions to the average 
{5T)n- The first of these arises from the diagonal ap- 
proximation for the average in the nth cavity. As before, 
we define the number m = ma = m/s as the number 
of times the trajectories a and f3 enter the nth cavity 
from the (n — l)st cavity. The first additional contribu- 
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FIG. 6: (Color online) Two examples of trajectory pairs con- 
tributing to {5T}n in the presence of time-reversal symmetry. 




tion to {6T)n then involves the pairing of segments a^+j 
with the time-reversed of l3u+v-j-i, i = 0,...,w — 1, and 
l<M<it-|-f<m, 



oo oo 



v — l m— u+l 



oo oo oo 



+^EE E 

u—1 v—l m—u-\-v-\-l 



m— 3 



(32) 



Examples for trajectory pairs contributing to the first 
and second line in Eq. (32) are shown in Fig. |6] The 
second contribution comes from small-angle encounters 
inside nth cavity involving the segments a„, au+v, 
and Pu+v 1 where l<u<u + v<m. For this contribu- 
tion one finds 



9c 



(33) 



are 



Examples of trajectory pairs contributing to {ST^^ , 

shown in Fig. [t] Using T = T — T"^ and adding these 
two contributions to Eq. (p2|) one finds 



{ST), 



-T{n - l)trr(n - 1) 



(34) 



up to corrections of order 1 / g1 . 

For the fluctuations of 5T in the presence of time- 
reversal symmetry one finds an extra contribution from 
pairing aj with the time-reversed of (3'^j^i_j, j = 1, . . . ,u, 
or a'j with the time-reversed of (3u+i-j, while following 
diagonal pairing rules for all other segments. Proceeding 
as before, we find 

{5T{p^,P2)5T{p'„p'^))^ = -[T{p',,p2)T{p,,p'^) 

9c 

+ 2gip,p[)g\p2,p'^)] 

+ 0(l/5?). (35) 



where 



G{p,p') = G{p',p) 

= [Si2Sl2S2l]{p,p'). 



(36) 



The stochastic process defined by Eqs. ( |34[ ) and (35 1 pre- 
cisely mirrors the stochastic process ^ for the quantum 
mechanical matrix T"^ = >S'^2*5'i2- 



FIG. 7: (Color online) Two examples of trajectory pairs con- 
tributing to {5T)„ in the presence of time-reversal symmetry. 
The trajectory pairs have a small-angle encounter in the nth 
cavity. 



C. Recursion relations for the moments of T 

The direct construction of recursion relations for the 
moments of T is an alternative to the stochastic approach 
that avoids extending the use of the properties (i) and 
(ii) of Sec. Ill to sample-specific quantities and the ne- 
cessity to define a statistical ensemble by varying the 
Fermi energy or shape of each cavity individually. The 
construction of recursion relations for moments of T pro- 
ceeds in the very same manner as the construction of the 
stochastic recursion relations for T, with the additional 
requirement that trajectories are piecewise paired in all n 
cavities, not only in the nth cavity. Since the arguments 
of the preceding section did not rely on the structure of 
the trajectories in the first n— 1 cavities, one immediately 
concludes that the recursion relations for the moments of 
T derived this way are identical to the recursion relations 
for the moments of T one obtains from the stochastic ap- 
proach. Starting from the stochastic recursion relations 
^ ~ ~ 351) 



(22) and (31 1 or (34) and ( |35[ ) one arrives at the same 
hierarchy of recursion relations ([7| derived for an array 
of disordered cavities. 

We illustrate this procedure for the recursion relation 
for the first moment (tr T(n)) in the absence of time- 
reversal symmetry. Following the rules of the semiclas- 
sical formalism, the average (tr T(n)) is determined by 
trajectory pairs a, f3 that are piecewise paired through- 
out the entire array of cavities. The trajectories can have 
small-angle self encounters, at which the pairing between 
a and (3 can be changed. Each pair of trajectories is clas- 
sified by the number m of times the trajectories enter the 
nth cavity from the (n — I)st cavity. The m segments of 
each trajectory in the nth cavity are labeled ai, . . . , 
and . . . ,/ 3m- 

As in Sec. IV A it will be sufficient to consider pairs 
of trajectories in which each segment aj is paired with 
the corresponding segment f3j, j — 1, . . . ,m. Trajectory 
pairs with self encounters in the nth cavity or with non- 
diagonal pairings in the nth cavity give contributions to 
(5(trT) that are a factor 1/gc smaller than the leading 
contribution. For the diagonal pairing, the remaining m 
segments of the trajectories a and P that reside in the 
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a = f) 



a = l3 






FIG. 8: (Color online) Examples of trajectory pairs con- 
tributing to the full average {T{n)) in the absence of time- 
reversal symmetry. The trajectory pairs have m = 1 (top), 
m = 2 (center), and m = 3 (bottom). The left panel shows 
trajectory pairs without small-angle self intersections; The 
right panel shows trajectory pairs with two small-angle self 
intersections in the first n — 1 cavities. 



first n—\ cavities precisely generate {tv'T(n — \)[ivTZ{n — 
1)]™"^), where the averaging brackets refer to variations 
of the Fermi energy for the entire array of cavities. Hence, 
we find 



(trr(n)) = £ 



(trT(n- l)[tr7^(n- 1)]™"^) 



ffctr T(n - 1) 
- tr7^(r^ - 1) 



(37) 



A schematic representation of trajectory pairs contribut- 
ing to Eq. ( 37 1 with up to two self-encounters in the first 
n — 1 cavities are shown in Fig. [8] Using unitarity to 
express TZ in terms of T and subtracting (trT(n — 1)), 
one then finds 



J(tr r(n)) (trr(n)) - (trr(n- 1)) 

= i((trr(n-l))2). 
5c 



(38) 



This is the same equation as one obtains from taking the 
average of the increment tr ST in the stochastic approach. 



V. FIELD THEORY FORMULATION 

In this section, we will approach the localization phe- 
nomenon from a different perspective. We will use that 



the quantum dot array depicted in Fig. [T] supports hier- 
archies of different types of field-fiuctuations in a field- 
theoretic description. These fluctuations reflect the fate 
of density distributions in classical phase space under the 
dynamical evolution of the system. Each of these fluctua- 
tions, thus, comes with a characteristic 'relaxation time', 
i.e., a time scale on which the fluctuation decays. (For ex- 
ample, fluctuations inhomogeneous in the sector of phase 
space describing an individual quantum dot will decay 
on a time scale comparable to the time of flight through 
the dot, etc.) In the description of low energy phenom- 
ena such as the zero frequency (DC) conductance, modes 
operating at short time scales can be treated perturba- 
tively. Their feedback into the sector of long time scales 
then stabilizes a 'low energy theory'. In the following, 
we will derive a theory that is minimal in that it con- 
tains information equivalent to that stored in the Fokker- 
Planck equation of localization. The strategy pursued 
here parallels one applied previously '^''J to the problem of 
dynamical localization in the quantum kicked rotor (also 
known as the "standard map" .) One difference is that the 
spectrum of different modes encountered in the present 
problem happens to be more complex. Our logics also 
resembles that of Ref. |52j where it had been shown that 
the relevant low energy theory of an ergodic quantum sys- 
tem contains the information otherwise stored in random 
matrix theory. 

Technically, our discussion will be based on a formu- 
lation of the array in terms of the ballistic nonlinear a- 
modcl.''^"^ A quadratic approximation in energetically 
high lying modes generates an effective low energy the- 
ory wherein each quantum dot is treated as a structure- 
less ('ergodic') entity. This theory will be equivalent to 
the celebrated nonlinear cr-model of disordered quantum 
wires,!^ a model that predicts Anderson localization on 
large length scales. We will see that the parameters sta- 
bilizing the hierarchical mode integration are the same 
as those utilized in previous sections of this paper. 

Conceptually, the hierarchical scheme is an alterna- 
tive to an indiscriminate perturbative integration over 
all modes in one go. That latter scheme would be essen- 
tially equivalent (see Ref.^^ for a discussion in the con- 
text of spectral statistics.) to a semiclassical expansion 
in terms of paired trajectories. In this sense, the hierar- 
chical mode integration processes the information stored 
in the statistics of trajectories by different means. 



A. Field theory of the quantum dot array 

Our starting point will be the description of the quan- 
tum dot array in terms of the supersymmetric ballistic 
nonlinear cr-model. This theory is obtained by averag- 
ing exact functional representations of Green functions 
over an energy interval of width AE centered around the 
uniform Fermi energy Ep of the array. A subsequent 
saddle point approximation (stabilized in the parameter 
Ep/AE ^ 1) then obtains a fleld theory in classical 
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phase defined by 



DT exp(-5[T]), 



(39) 



with 



S[T] = ^ |^(dx)tr(r*A{i/,r-i} 



^rcg[T]. 



(40) 



The integration variables in (|39|, T(x) = {T"" (x)} 
are (super) matrix valued fields defined on shells F = 
{x|iJ(x) — Eq} of constant energy in classical phase 
space. Further, x = (q, p) where q and p are coordi- 
nates and momenta, respectively, H{yi) is the Hamilto- 
nian function of the system, the integral over the energy- 
shell is normalized to the (spatial) volume of the sys- 
tem, Jp((ix) = Vol, and v is the single particle den- 
sity of states per volume, v = l/(AVol). For time 
reversal and spin rotation invariant systems (orthogo- 
nal symmetry class, f3 — I), the 'internal' structure of 
the matrices T"" is described by a composite index 
a = {a,r,t), where a = +/— discriminates between the 
advanced and retarded sector of the theory, r = b, f dis- 
criminates between commuting and anticommuting sec- 
tors, and t — 1,2 accounts for the operation of time re- 
versal. Time reversal symmetry reflects in the relation 
(TT'^r~^)(q, — p) = T~^(q, p), where r is a fixed matrix 
whose detailed structure will not be of concern through- 
out. For time reversal non-invariant systems (unitary 
symmetry, /3 = 2) no time-reversal structure exists and 
a — (a, r). In either case, the matrices T carry a coset 
space structure in the sense that configurations T and 
TK are to be identified if [K,A] = 0, where A = erf 
and "ar" stands for action in advanced/retarded space. 
Finally, the regulatory action 



'5'rog [T] 



; J (dx) str^AT-i * at), 



5i0 



determines the 'causality' of field fluctuations, but will 
otherwise not be of much relevance throughout. 



The fluctuation behavior of the fields T in ( 39 1 is gov- 
erned by the classical Liouville operator {H^} (where 
{ , } is the Poisson bracket.) Quantum mechanics en- 
ters the pro blem through the presence of Moyal products 



in (39 1. In essence, this product operation limits 



the maximum resolution of the theory in phase space to 
scales of the order of a Planck cell. In the following sec- 
tions we reduce the above 'bare' theory to an effective 
low energy theory describing localization phenomena. 



B. Hierarchical mode integration 

Before turning to the technicalities of the mode inte- 
gration, let us describe the relevant hierarchies in quali- 
tative terms: fluctuations inhomogeneous in the phase 
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FIG. 9: (Color online) Top left: cartoon of an array of weakly 
coupled quantum dots. Center left: cut region containing two 
dots. The shaded region indicates the real space support of 
the coupling phase space Tc- Bottom left: plot indicating the 
profile of weakly fluctuating fleld configurations: constancy 
throughout each dot, gradual variation in the coupling re- 
gion. Center right: zoom into the coupling region. The phase 
space region Tc contains points x = (q, n) that will migrate 
between the two dots in short time, i.e., without undergoing 
back scattering. (The opening angles represent those direc- 
tions n that meet this criterion.) Bottom right: schematic 
plot of real space profile of phase space fiuctuations generat- 
ing the inter-dot coupling. 



space sector representing individual dots are expected 
to relax on short time scales comparable to the time 
of flight across the dot, On the other hand, rel- 

ative fluctuations in the configuration of different dots 
can survive up to time scales of the order of the dwell 
time ^ tf. To describe this hierarchical decay profile 
within our field theory framework, we focus on a sec- 
tion of the array containing two neighboring quantum 
dots (cf. Fig. [9|. The fields T(x) representing phase 
space fluctuations in this subsystem may be decomposed 
as r = TgTf, where Tsj are 'slow' and 'fast' fluctuations, 
respectively. The slow fluctuations are (i) homogeneous 
within each dot separately. In particular, (ii) they do not 
vary in any spatial cross section transverse to the array, 
and (iii) do not depend on momentum. However, (iv) 
the weakness of the inter-dot coupling (td 2> t/) leaves 
room for gradual fluctuations of the slow modes as we 
pass from one dot into the other. This suggests a pa- 
rameterization T's(x) = Ts(q), where q is the component 
of q parallel to the longitudinal direction of the two-dot 
system. Point (i) above means that Ts{q <^ 0) = and 
Ts{q ^ 0) = Tr, where Tl^r are the constant slow mode 
configurations of the left and the right dot, respectively. 
As we are going to check in a self consistent manner, (v) 
relative fluctuations between left and right dot are sup- 
pressed (in a parameter of the order of the number of 
transverse quantum channels supported by the connec- 
tor region) so that a leading order expansion in relative 
fluctuations Tj^T^^ is sufficient. (Conceptually, this ex- 
pansion is equivalent to the Kramers-Moyal expansion 
employed in the derivation of the Fokker-Planck equa- 
tion above.) 
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The field Tf encapsulates all other fluctuations, i.e., 
fluctuations that do not meet the criteria (i-v). Generic 
fluctuations of this type - think of fast fluctuations deep 
inside the phase space of one of the dots - are strongly 
gapped and do not couple to the slow fluctuations. (For- 
mally, this decoupling manifests itself in an effective 'or- 
thogonality' between the modes representing these fluc- 
tuations.) However, there is one sector of phase space, 
Tc, in which fast and slow fluctuations talk to each other. 
The domain Tc includes all points in phase space which 
pass from one dot to the other in a time of order tf <C td , 
much smaller than the typical dwell time (cf. the shaded 
area in Fig. [9]). This region is special in that it overlaps 
with the domain of gradual variation of the slow fields 
(cf. point (iv) above.) As we will see, a perturbative in- 
tegration over fast fields in F^ effectively determines the 
slow field coupling between the two dots. 

To prepare the integration over the fast fields in Tc, 
we need to bring the notation to a more explicit level: 
assuming that the bulk dynamics is ballistic, {H, } — 
vp'n. • V where n is the unit vector in momentum space, 
and vp the Fermi velocity. Current conservation in the 
specular refiection at the system boundaries translates to 
the effective boundary condition /(x, n) = /(x, n) , where 
/ is a phase function subject to the action of {H, }, x a 
boundary point and n the direction vector with flipped 
normal component. 

We next parameterize fluctuations as Tsj — 
ex.p{Wsj), where the field generators carry a block struc- 
ture (in advanced/retarded space). 



(For further technical details on this representation we 
refer to Ref.1551) We now substitute these generators into 
the action, expand to leading order in Wf, and integrate. 
This leads to the effective action 

SaSm = I (^[s^'HTs,Wf]y'j + 5,eg[T,], (41) 

where (...) = J DWf exp{-S'^^'>[l,Wf]){. . .), and 
^Ts, Wf] is of nth order in Wf. Due to the isotropy of 
Ts in momentum space and the linearity of the Poisson 
bracket in n, the action does not contain a contribution 
of zeroth oder in the fast fields. The dominant contribu- 
tion to the coupling between fast and slow fluctuations 
is given by the linear term, 

S^^'>[T,,Wf]^ f3TThi^VF I dqdnn||Str[W^,(q,n)$(g)] 

(42) 

where we have introduced the abbreviation $ = (dqTs) * 
ATjT^, the integration over the direction of momentum is 
normalized as J dn = 1, riy is the component of n parallel 
to the longitudinal direction and the Moyal product has 
been omitted (which is permissible due to the general re- 
lation J {dx){f * g){x.) = /((ix)(/ (7)(x) for the integrated 
product of two functions - presently, matrix elements of 
<& and W - in phase space.) 

Neglecting contributions of 0{^Wf) (as compared to 
the C'(<I>W/)-terms above) the quadratic Ws-action is 
given by 



S^^^[T,,Wf] ~ S^^^[l,Wf] 

= [jirhv i dq(instr(_B/(q, n)(ui;-n • V + (5)i3/(q, n)). 

Fast field fluctuations can now be integrated out according to the prescription 
(str(A(q, n)B/(q, n)) str(^(q', n')5/(q', n'))) = <5(n - n')n(q, q'; n), where 

n(q,q'; n) = -^{vpn ■ V + 5)-\ci, q'). 

Specifically, the effective action is given by 

5cff[r,] = {p^huvpf I dqdq'dnnfstr ($2i(g)*i2(g))n(q,q';n) 

So far, we have not made reference to the specific properties of the phase space region Tc. We now assume that the 
corridor connecting the dots has wave-guide properties, in that it (a) does not contain significant backscattering, and 
(b) the restriction of the Liouville operator to Tc has plane-wave like eigenfunctions characterized by a conserved 
longitudinal/ trans verse momentum \i±/k. We also assume that (c) the slow fields smoothly interpolate between 
and Tfj in a region centered around the longitudinal coordinate q — Q. The presumed proximity — Tr implies that 
<i>((7) ~ dqW{q)K can be linearized, where we suppressed the slow field index 's' for notational transparency. Under 
these circumstances, and noting that the integration over qj^ implies a projection onto the zero-momentum sector 
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kj^ = 0, we obtain 



dk e*'=(9-«') 



J-c J J 2mvFn\\k + d 

^-PttHvScVf J dqdq'str{^2i{q)^i2{q')) J dn\n\\\e ' ^ ^ 



— const. X hi/ScVFStT{{BR ~ Bl){Bb. — -Bl)), 



where the indices refer to "ar" -space, const. — 0(1) is a 
constant, Sc the transverse cross section of the contact, 
and i?ij,L are the generators of the slow fields in the left 
and the right dot, respectively. Noting that the density 
of states per volume, v ^ mkp~^ (where kp — vpin is 
the Fermi momentum) and Sck'j~^ N is proportional 
to the number of transverse channels, TV, supported by 
the connector region, the prefactor can be written as 
const. X iV 3> 1, a number which we assume large lest 
a semiclassical description of the contact becomes mean- 
ingless. We also note that the quadratic form may be re- 
placed by its unique rotationally invariant generalization 
to the full field manifold, str((i?fl^ - Bl){Br ~ Bl)) 
\siv{QLQFi), where Q — TAT^^. While a quadratic ex- 
pansion of the latter expression reproduces the bilinear 
term, the largeness of N implies that typical values con- 
tributing to the ^-integration, B ~ N~^^^ <C 1, which 
means that non-linear contributions become inessential 
in the limit of large channel numbers. We thus conclude 
that the coupling term can be rewritten as 

SeslQ] = const. X Nstr{QLQR)- 

Finally, the obvious generalization of the above two-dot 
construction to an entire quantum dot array reads as 



ScsiQ] = const. X ^str(Q,„(5„+i), 



• Semiclassically speaking, a 'mode' represents the 
coherent propagation of a retarded and an ad- 
vanced Feynman amplitude along classical trajec- 
tories in phase space. Locally, the semiclassical dy- 
namics of such composites is described by the Li- 
ouville operator, as is manifest in the action (39 1. 



This trajectory interpretation helps in identifying 
the relevant fast modes. E.g., in the system de- 
picted in Fig. [9] the coupling between the ergodic 
slow modes of each dot, Qmi is dominated by tra- 
jectories swiftly propagating from one dot to the 
other, i.e., modes emanating at phase space points 

xeFe. 



• Although this identification of fast modes rests on 
specific model assumptions (no backscattering in 
the contact region, etc.), the result (43) is rea- 
sonably universal. For example, a somewhat more 



(43) 



elaborate construction will show that the same ac- 
tion describes connector regions containing chaotic 
scattering and momentum relaxation by them- 
selves. (The latter complication would manifest 
itself in an altered value of td, though.) Gener- 
ally speaking, the low energy physics of the system 
will reduce to ScS, as long as the dots are isolated 
from each other in the sense tf <^ td. 



where Qm is the Q-matrix representing the m-th dot and 
we used that the channel number N ^ (tdA)~^. Before 
turning to the discussion of localization properties, a few 
remarks on the construction above are in order. 

• Conceptually, the above reduction programs in- 
volves three steps: 1) identification of 'low energy 
modes', i.e., modes that decay on the largest time 
scales of the problem, 2) identification of 'high en- 
ergy', or quickly decaying modes, and 3) pertur- 
bative integration over those fast modes that will 
conceivably couple to the slow modes. In principle, 
that integration can be explicated for any fluctua- 
tion in the problem. In practice, however, only few 
modes will effectively couple to the slow degrees of 
freedom, and these relevant fluctuations are best 
identifled by semiclassical considerations: 



C. Localization from the effective action (43) 



The effective action (43) is equivalent to a 'lattice ver- 



sion' of the diffusive nonlinear cr-model of quasi-one di- 
mensional disordered wires. Indeed, we may pass to a 
continuum limit 



tdA 



tdA 



dx sir{dQdQ), 



where Qm Q{x) is replaced by a smooth field, x — ma, 
and a the spacing between the dots. Comparing to the 
standard form of the diffusive model,'^ where the action 
is ~ ^ / dxstridQdQ) with ^ the localization length, we 
are led to the identification — ^ ~ 
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VI. CONCLUSION 

In the preceding sections we showed how a theory of 
Anderson locaHzation can be constructed for a sample 
in which the microscopic electron dynamics is ballistic, 
rather than disordered-difFractive. Our theory of "ballis- 
tic Anderson localization" paraphrases the scaling ap- 
proach to localization in disordered qu antum wires of 
Dorokhov, Mello, Pereyra, and Kumar j^S129| using the 
language of the trajectory-based semiclassical formalism. 
As noted in the introduction, the interest of construct- 
ing such a semiclassical theory is not the structure of the 
theory itself or the phenomena it explains. Like most 
semiclassical theories of quantum corrections in the per- 
turbative regime, the structure of the theory closely re- 
sembles the structure of its fully quantum mechanical 
counterpart for disordered metals, whereas the observed 
phenomena are the same in the ballistic and disordered 
cases. Instead, the interest of the semiclassical theory 
is that it shows how quantum effects that were known 
from disordered metals arise if the electron dynamics is 
ballistic. 

On hindsight it should not come as a surprise that 
a theory of ballistic Anderson localization can be con- 
structed by adapting the derivation of the Dorokhov- 
Mello-Pereyra-Kumar equation for a disordered quan- 
tum wire. This point is best made by reconsidering 
Dorokhov's original derivation.^^S! In this derivation, im- 
purity scattering is treated in the Born approximation. 
All quantum mechanical amplitudes are squared into 
quantum probabilities. Hence, it is sufficient if one can 
replace the quantum mechanical probabilities by classical 
ones. Such a replacement is a standard procedure when 
connecting quantum mechanical and semiclassical theo- 
ries. Its implementation for the array of chaotic cavities 
is what is done here. 

It is interesting to observe that, while small-angle en- 
counters form a crucial link in our understanding of quan- 
tum interference corrections in ballistic conductors, their 
role is very limited in our description of localization in 
quasi-one dimension: They serve to cancel spurious con- 
tributions from trajectories that enter the last cavity 
of the array more than twice. (In Dorokhov's original 
approach such processes are excluded automatically be- 
cause of the condition that the length of the wire is in- 
creased by an am ount SL much smaller than the mean 
free path iW^^ Implicitly, small-angle encounters do 



play a much more important role in our theory, however, 
because they help to preserve unitarity in the semiclas- 
sical theory. Unitarity is a key ingredient of both the 
quantum-mechanical derivation of the DMPK equation 
and the present semiclassical derivation. We note that 
unitarity has played an important role in other exten- 
sions of the semiclassical framework beyond its previously 
assumed domain of validity: it is used to relate weak 
localization and enhanced backscattering, thus enabling 
a semiclassical description of wea k localiz ation without 
reference to small- angle encounters ,12^1211111 and it is used 
to obtain an alternative expression for the spectral form 
factor, allowing its calculation in the non-perturbative 
regime by considering periodic orbits of duration below 
the Heisenberg time only.*^ 

We have also shown that the dynamical information 
entering the semiclassical approach can be processed by 
different means to derive an effective low energy field the- 
ory of the system. This latter approach is based on the 
concept of 'modes', fluctuations in phase space decaying 
on parametrically different time scales. A successive in- 
tegration over short lived modes stabilizes an effective 
action of the most persistent modes in the fluctuation 
spectrum. In the present context, that low energy limit 
turned out to be equivalent to the diffusive nonlinear a- 
model of disordered wires. While the field theory ap- 
proach is arguably less explicit than the direct classifica- 
tion of trajectories, it enjoys the advantage of high com- 
putational efficiency, a paradigm previously exemplified 
on the problem of universal spectral correlations."''^ For 
example, the above mentioned condition of unitarity, as 
well as the symmetries relating trajectories to their time 
reversed are built into the field theory approach from 
the outset; there is no need for explicit bookkeeping in 
terms of encounter processes. The price to be payed for 
this compactness in the description is a higher level of 
abstraction, though. 
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